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SUMMARY 

The trajectories of droplets in the air flowing past a 36.5-percent- 
thick Joukowski airfoil at zero angle of attack were determined. The 
amount of water in droplet form impinging on the airfoil, the area of 
droplet impingement, and the rate of droplet impingement per unit area on 
the airfoil surface were calculated from the trajectories and cover a 
large range of flight and atmospheric conditions. With the detailed 
impingement information available, the 36.5-percent-thick Joukowski air- 
foil can serve the dual purpose of use as the principal element in instru- 
ments for making measurements in clouds and of a basic shape for estimating 
impingement on a thick streamlined body. 

Methods and examples are presented for illustrating some limitations 
when the airfoil is used as the principal element in the dye-tracer 
technique. 


INTRODUCTION 

A commonly used technique for making measurements in icing clouds, 
either in flight or in icing tunnels, involves the collection of ice on 
cylinders of different diameters exposed to the airstream. This technique 
is known as the rotating multicylinder method and is described in refer- 
ence 1. Another technique for measuring cloud droplet sizes principally 
in tunnels and with temperatures above freezing is described in reference 
2 as the dye -tracer technique. In this method, a blotter-wrapped body, 
such as a cylinder or airfoil for which theoretical droplet impingement 
data are available, is exposed to an airstream containing a dyed-water 
spray cloud. The amount and local distribution of dye residue in the 
blotter are analyzed and are compared with theoretical values in order to 
obtain the sizes of the droplets in the impinging cloud. 
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Several other methods involving photography, sampling techniques, 
light scattering, droplet impaction, or combinations of two or more of 
these principles, such as are described in reference 3 (see also refer- 
ences listed in ref. 3), have been employed in the process of developing 
cloud measuring instruments. The scope of this report is confined to 
presenting data and discussion intended to refine the use of the dye- 
tracer technique. 

The impingement data obtained from droplet trajectory calculations 
are required for either the rotating multicylinder method or the dye- 
tracer technique. The calculation of droplet trajectories requires a 
knowledge of the airflow field surrounding the aerodynamic body being 
studied. In order that the trajectory calculations apply in principle, 
the calculated air velocity components surrounding the body must be the 
same as those existing in the atmosphere during the measurements. The 
velocity components used in the calculations are usually obtained from 
the solution of theoretically derived flow equations for a body immersed 
in an ideal incompressible fluid. The equations do not account for 
viscosity, separation of the flow from the body, effects of compressibility, 
turbulence, and other factors. The effects of compressibility of the air 
on droplet trajectories have been found to be negligible up to the flight 
critical Mach number of the body (ref. l). With most well designed aero- 
dynamic bodies and when used for measuring atmospheric cloud droplets of 
interest in icing, the effect of viscosity is small as it is confined to 
the boundary- layer region. 

A concern for the effect of flow separation around cylinders is 
expressed in reference 2, in which a comparison of the theoretical and 
wind tunnel experimental surface velocities on cylinders is given. The 
comparison is repeated herein as figure 1. It is apparent that the ratios 
of local surface velocity to free-stream velocity differ over the entire 
cylinder surface. The difference between experiment and theory may be 
due to early separation of the flow. This separation may induce a flow 
pattern ahead of the cylinder that is not attendant to a right circular 
cylinder but rather to some virtual shape related to the position of 
separation on the cylinder surface. The large differences in the local 
velocity on the cylinder surface indicate that the airflow field ahead 
of the cylinder may be considerably different from the theoretical flow 
field. Consequently, actual droplet trajectories may be different from 
the theoretical calculations of reference 1 and may contribute to errors 
and to such anomalies as discrepancies in the cloud measurements (droplet- 
size distribution and liquid-water content) when the size of the cylinder 
is varied. These discrepancies are related to cylinder size in such a 
manner as to be inconsistent with the theoretical analyses of reference 
1. The varied deviations of measured surface velocity among different- 
sized cylinders from the theoretical surface velocity (fig. l) imply an 
explanation for the cylinder-size trend. 
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The experimental work of reference 2 suggests the use of a "blunt- 
nosed body such as the cylinder but with a streamlined trailing section 
that minimizes flow separation. The blunt-nosed shapes are preferred, 
especially when used with the dye-tracer technique described in reference 
2. The principal reasons for this preference are that the limits of im- 
pingement and distribution of water on the surface are more easily measured 
on bluff bodies . 

On the basis of wind tunnel experience, a thick Joukowski airfoil 
section was investigated as a possible streamlined, blunt-nosed body. A 
Joukowski airfoil has many of the desirable geometric qualities and, in 
addition, is amenable to simple calculations of the flow field. A 60- 
percent-thick Jouskowski airfoil section was tested in the wind tunnel 
and was found to have serious airflow separation. A section arbitrarily 
reduced to a thickness of 36.5 percent was a satisfactory bluff body with 
only a small difference between the theoretical and experimental surface 
velocities (fig. 2). This airfoil gives a surface velocity ratio within 
4 percent of the calculated ideal. When the air velocity distribution on 
the 36.5-percent-thick airfoil was found satisfactory and the airfoil 
retained the desired bluntness, no other airfoil was examined. 

Droplet trajectories were calculated for the 36 .5 -percent -thick 
Joukowski airfoil at the NACA Lewis laboratory. The calculations help 
explain some of the difficulties, which are discussed later, encountered 
in the use of cylinders. With detailed information on droplet impingement 
available, the 36 .5 -percent -thick Joukowski airfoil serves a dual purpose. 
Besides providing an aerodynamic body for replacing cylinders when they 
exhibit undesirable surface velocity distributions in cloud measuring 
instruments, the calculations serve as an immediate contribution of in- 
formation on impingement on very thick airfoil sections used in some 
applications such as radar housings. Streamlined struts closely resembling 
thick-sectioned Joukowski airfoils are also often used for mounting in- 
struments, fuel tanks, and armament projecting into the airstream. 


CALCULATIONS OF DROPLET TRAJECTORIES 

In order to find the rate and distribution of droplet impingement 
on the surface of a body, it is necessary to determine the cloud-droplet 
trajectories with respect to the body. The method used for the 36.5- 
percent-thick Joukowski airfoil in calculating the droplet trajectories 
is essentially the same as described in reference 1 for the cylinder. 
Assumptions necessary for the solution of the problem are: (l) At a large 

distance ahead of the body the droplets are at rest with respect to the 
air; (2) the only external force acting on the droplets, as the airfoil 
approaches the droplets, is the drag force due to the relative velocity 
of the air with respect to the droplets; and (3) the droplets are always 
spherical and do not change in size. 
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Differential Equations of Droplet Motion 

The differential equations that describe the motion of the droplets 
are obtained by equating the drag force with the rate of change of momentum 
of the droplet. The equations apply to the motion of droplets in a plane 
coincident with the airfoil section as shown in figure 3. The airfoil 
location in the coordinate system used in other MCA reports (such as ref. 
4) is retained herein. The geometric chord line of the airfoil is col- 
linear with the x-axis of the rectangular coordinate system, and the lead- 
ing edge is placed at the origin of the coordinates. At an infinite dis- 
tance ahead, of the airfoil, the uniform airflow carrying the cloud drop- 
lets is assumed to be approaching the airfoil from the negative x-direction 
and parallel to the x-axis . 

The equations of motion, as derived in reference 1, are 

dv x C lP e 1 , , 

dT = 24 K ' U x “ V 

( 1 ) 

^ i T ) 

dT 24 K ' U y V y' 

where Re is the. local Reynolds number with respect to the droplet diam- 
eter 2a and the local relative velocity between the air and the droplet. 
(All symbols are defined in appendix A.) Thus, 



The coefficient in equation (2) is called the free -stream Reynolds number 
ReQ . Hence, 



The dimensionless number K, the inertia parameter, is defined as 


The coefficient of drag Cp is obtained from experimental data as a 
function of the local droplet Reynolds number Re (ref . 5) . 


The equations are expressed in dimensionless form in order to main- 
tain the number of calculations at a minimum and to simplify the presenta- 
tion of the results. The dimensionless coordinates x and y are ratios 
of the actual distance to the chord length L, which is the unit of 
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length. The dimensionless air velocity components u x and Uy and the 
droplet velocity components v x and Vy are ratios of the actual veloc- 
ity to the free-stream velocity, which is the unit of velocity. Time is 
expressed in terms of the dimensionless quantity t = tV/L. The unit of 
time is the time required to go a distance L at a speed V. 

An examination of equations (l) and (2) shows that the characteris- 
tics of the trajectories depend only on the values of the two dimension- 
less parameters K and ReQ. Thus, a unique droplet trajectory exists 
for each combination of K and Re 0 . 


Method of Solution 

The differential equations of motion (eqs. (l) ) are difficult to 
solve because values of the velocity components and the factor containing 
the coefficient of drag depend on the position and velocity of the droplet 
at each instant and, therefore , are not known until the trajectory is 
traced. The values of these quantities must be fed into the equations as 
a trajectory is developed. This was accomplished by using a mechanical 
differential analyzer constructed at the HACA Lewis laboratory for this 
purpose (ref. 6). The results were obtained ih the form of plots of 
droplet trajectories in the coordinates of figure 3. 


The velocity field around the airfoil, required for the solution of 
the trajectory equations, was obtained from equations for incompressible, 
ideal flow around a Joukowski airfoil such as given in reference 7. A 
simple translation of coordinates and scale is necessary when the usual 
equations for the Joukowski airfoil (ref, 7) are used in order to place 
the airfoil as shown in figure 4 with one chord as the unit of length. 

The surface coordinate points are presented in table I, and the airfoil 
shape is illustrated in figure 4. The relation between chordwise and 
surfacewise positions, which is often very useful, is given in figure 5. 

The equations of motion were solved for Various values of the param- 
eter K. The inertia parameter K is a measure of the droplet size, the 
flight speed and size of the airfoil, and the viscosity of the air through 
the relation 

K = 1.704X10" 12 d 2 U/|iL (5) 

The density of water, 1.94 slugs per cubic foot, is included in the con- 
stant. For each value of K, a series of trajectories was computed for 
each of several values of free-stream Reynolds number 

Re 0 = 4 . 813x10 ~ 6 dp a U/p 


( 6 ) 
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A convenient graphical means for determining values of the dimensionless 
parameters K and ReQ in terms of airplane speed (free-stream veloc- 
ity), chord length, altitude, and droplet size is presented in appendix 
B of reference 4. 


A dimensionless parameter <P 

_Re§ 
<P = "X" 


defined as 


= 13.594 



(?) 


is also used in presenting the data. The parameter <P is valuable in 
that it is not a function of droplet size. The parameter <P is an im- 
portant concept in the interpretation of cloud measurements in which the 
droplet size is not measured directly and is an unknown that must be cal- 
culated (see ref. l) , The use of any two of the three dimensionless 
parameters K, ReQ, or <p will completely define a data point with re- 
spect to the flight and meteorological conditions. 


IMPINGEMENT RESULTS 


For a symmetrical airfoil at zero angle of attack, the airflow field 
for negative values of the y- ordinate is a mirror image of the flow field 
around the upper airfoil surface. Thus, only the impingement on the upper 
surface is presented. 


The rate of water impingement and the manner in which it is distrib- 
uted on the surface of the airfoil can be obtained if the starting point 
of a droplet trajectory is known with respect to the point of impingement 
on the surface. The starting ordinate yQ at infinity of any impinging 
trajectory can be found in figure 6 with respect to the point of impinge- 
ment on the surface. The values for the starting and ending positions of 
the trajectories are shown in figure 6 for five values of free-stream 
Reynolds number. For each value of ReQ, curves for several values of 
K are given. Since the curves for very low values of K are difficult 
to read from the scale used in presenting figure 6, the same data are 
plotted in figure 7 for three low values of K. 

The calculated points for the curves are shown in figures 6 and 7 in 
order to give an appreciation of the precision of the curves. Informa- 
tion on precision is often necessary when the airfoil is used as the 
principal element of an instrument. Because the mechanical analogue used 
for the calculations is substantially a graphical method of solution, the 
data points in figure 6 have a tolerance of ±0.001 on both values of s 
and values of yQ. In figure 7 the tolerance is ±3.0005 for s and 
±0.0002 for yQ. The difference in tolerances is caused by the change in 
scale factors during the calculations. The tolerance for the values of 
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s along the dashed line in figure 6 (s may ) described in following sec- 
tion) is ±0.002 because the determination of the point of tangency pro- 
vided an added factor of uncertainty. 


Extent of Impingement 

The limit of impingement is determined by the point of tangency on 
the airfoil surface of the tangent trajectory. In figure 6 the limits of 
impingement s max lie on the dashed lines, which are the loci of the 
termini of the curves and which determine the point of impingement of the 
tangent trajectory for each combination of K and EeQ. For convenience, 
the rearward limits are replotted in figure 8 as a function of K for 
various values of ReQ. The distances are measured on the airfoil sur- 
face from the point of intersection of the geometric chord line with the 
leading edge (fig. 3) in terms of the chord length. When the airfoil is 
used as the principal element of an instrument, the values of s max as 
a function of K for various values of <P are often more useful. This 
is shown in figure 9. 


Rate of Total Water Interception 

The rate of total water interception, in pounds per hour per foot of 
wing span, is determined by the spacing of the tangent trajectories (fig. 
3), by the speed of the aircraft, and by the liquid-water content in the 
cloud. The flight speed and size of the airfoil, as well as the droplet 
size in the cloud, are the principal variables that affect the spacing 
between the upper- and lower-surface tangent trajectories. For a sym- 
metrical airfoil at zero angle of attack, the water that strikes the air- 
foil is divided equally between the top and bottom halves and, therefore, 
is proportional to twice the spacing between the x-axis and the ordinate 
y 0 -|- an at infinity of the upper surface tangent trajectory. In figure 

6, the values of yQ ^. &n and those for s lie on the dashed lines. 
The rate of total waiter interception per unit span of the airfoil on that 
portion of the airfoil surface bounded by the upper and lower tangent 
trajectories can be calculated from the relation 

“max , 0.66U„Ly 0;tan (8) 

The rate of total water interception per unit of span can also be 
obtained from the collection efficiency E, which is defined as the ratio 
of the water contained in the droplets intercepted by the airfoil to the 
total water in the volume swept out of its path by the airfoil. On the 
basis of collection efficiency E = y^ tan /0.182, 

w max = O.lSUwLE 


(9) 
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The values of E are given in figure 10 as a function of K for various 
values of ReQ. Again, because of the usefulness in the application, the 
collection efficiency is repeated in figure 11 for various values of <p . 


Rate of Local Droplet Impingement 

The local impingement rate on the airfoil surface is often desirable 
knowledge in the design of thermal anti-icing systems or when the airfoil 
is used as a principal element of an instrument. The rate of water im- 
pinging between any two given points on the airfoil surface may be found 
by applying the results of figures 6 and 7 in the relation 

W = 0.33UwL(y 0A - y^ 2 ) (10) 


The local rate of droplet impingement per unit area of airfoil surface 
can be determined from the expression 


W c 


0.33UW 


% 

ds 


0 . 33Uw(3 


( 11 ) 


The values of the local impingement efficiency (3 as a function of the 
airfoil distance s are given in figures 12 and 13. These values were 
obtained from the slopes of the curves in figures 6 and 7, respectively. 

The values of p given in figures 12 and 13 are estimated to be in 
error by less than ±1 percent. Since the rate of total water impinging 
is directly related to the area under the P curves, a check on computa- 
tional accuracy of the values of p in figures 12 and 13 was also ob- 
tained by comparing the area under each p curve with values of collec- 
tion efficiency in figure 10. The area values checked within ±0.3 percent 
of the corresponding values for the rate of total water interception. 


DISCUSSION OP RESULTS 


Comparison of Impingement on 36.5- with 15 -Percent -Thick 
Joukowski Airfoils and Cylinder 

A comparison of available data on collection efficiency and limit of 
impingement for three bodies that vary in thickness ratio is helpful when 
interpolation as a function of thickness ratio is necessary. Data for a 
15 -percent -thick Joukowski airfoil are presented in reference 8 and for 
the cylinder in reference 1. The cylinder is included for comparison 
because it can be considered as a symmetrical airfoil with a thickness 
ratio of 100 percent. 
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The collection efficiency for the 36 .5 -percent -thick Joukowski air- 
foil is presented in figures 10 and 11 in terras of K, ReQ, and <p. 
Langmuir has suggested the use of a modified K parameter (refs . 5 and 
9) . The advantage of this modified parameter is that the data in terms 
of K and Reg or <p are essentially reduced to a single curve. For 
any given set of operating and meteorological conditions, which are rep- 
resented by one value of the modified parameter, the impingement charac- 
teristics of several bodies of different shape can be easily compared. 

The modified K parameter is known as the Kq parameter (see ref. 
9) which is defined as 


where X is the range a droplet would have when projected into still air 
if the drag is based on experimentally determined values (also known as 
true range) , and X g is the range of a droplet obeying Stokes 1 law of 
air drag for a sphere. The total collection efficiency E is given as 
a function of Kq in figure 14 for 15- and 36.5-percent-thick Joukowski 
airfoils and a cylinder. The values of E as' functions of Kq are shown 
in figure 14 as a single-valued curve for each body. Although no proof of 
the significance of Kq is available at this time, careful calculations 
reduce a family of curves, such as collection efficiency as a function of 
K and ReQ shown in figure 10, to a very narrow band of curves which 
essentially can be presented as a single-valued curve for purposes of 
comparison. For any value of Kq, which represents the same operating 
and meteorological conditions for all three bodies, the collection effi- 
ciency decreases with increase in thickness ratio. 

A comparison of the limit of impingement s„„„ is shown in figure 
15. The extent decreases with decrease in thickness ratio except for low 
values of Kq . Low values of Kq represent conditions of small droplets 
combined with large chord lengths and low speed. As an example for orien- 
tation, a value of Kq of 0.078 could represent an altitude of 10,000 
feet, a chord length of 18 inches, a speed of 200 miles per hour, and a 
droplet diameter of 18 microns . 


Impingement in Clouds of Nonuniform Droplet Size 

The data presented in figures 6 to 13 apply directly only to flights 
in clouds composed of droplets that are all uniform in size. The droplets 
in a cloud, however, may have a range of sizes. The design of equipment 
for protection against ice formation can be in error when based on a 
single droplet size such as the volume-median if a distribution of droplet 
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sizes is present in the cloud. The limit of impingement on the surface 
is governed principally hy the large droplets. Also, the total rate of 
impingement can vary considerably with changes in the size distribution 
of the cloud even though the liquid-water content and volume-median diam- 
eter are maintained. When the airfoil section is used as an element of 
a measuring instrument, the data in figures 6 to 13 must always be applied 
in a modified (weighted) manner. The method of weighting the data is 
discussed in detail in reference 2 in the application of cylinder data to 
the dye-tracer technique and is reviewed briefly, along with the dye- 
tracer technique, in appendix B herein. 


USE OF AIRFOIL AS SENSING ELEMENT 

An evaluation of the 36 .5-percent-thick Joukowski airfoil as a sensing 
element in an instrument can be made with illustrating examples . Much of 
the discussion on the airfoil also applies to other geometric shapes, par- 
ticularly cylinders, except that the range of limits of K and ReQ may 
vary. 


Basic Requirements of a Sensing Element 


A necessary requirement for a sensing element is that the unknown 
being measured by the element must have a well defined variation with 
respect to a known independent quantity. For example, the droplet size, 
which is the unknown being measured by the joukowski element, must be a 
well defined function of either surface extent or collection rate (either 
total or local, depending on which is used) , which presumably is obtained 
by measurement. For example, according to the curves in figures 9 and 
11, a well defined function exists between the limit of impingement (fig. 

9) or collection efficiency (fig. ll) and the inertia parameter K, which 
contains the factor of droplet size, for values of K from 0.04 to 5. 
Flight and meteorological conditions represented by values of K smaller 
and larger than the limits stated reduce the sensitivity of the thick 
Joukowski airfoil as a sensing element. Unfortunately, values of K less 
than 0.04 are often encountered in icing tunnel and flight operations. 

For example, for operation at 200 miles per hour, an altitude of 10,000 
feet, and an airfoil chord length of 18 inches, K is less than 0.04 when 
droplet diameters are less than about 8 microns. The upper limit, that 
is, K greater than 5, represents droplets greater than 87 microns provided 
the other conditions do not change. 

The preceding example illustrates the desirability of using several 
airfoils with different chord lengths (see ref. 10, appendix B) in measur- 
ing clouds with wide distributions in droplet size. A change in airfoil 
chord length has a large effect on the limits of sensitivity because the 
chord length L appears in the expression for the inertia parameter K 
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as well as in the determination of other factors' such as surface extent. 

If a 6-inch airfoil had been used in the preceding example, the values 
of K would have been three times larger (values of Reg are not affected 
by chord length) , and the readings for either surface distance of impinge- 
ment (figs. 8 and 9) or collection efficiency (figs. 10 and ll) would have 
been shifted into a region more favorable for the smallest droplets. This 
shift, however, hampers the large droplet measurements for which large 
airfoils are required. 


> 

> 
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Sensitivity for Droplet-Size Discrimination 

In the dye-tracer technique and other methods, the determination of 
the droplet-size distribution requires a known correlation between droplet 
size and limit of impingement. For sensitivity, a measurable variation 
between droplet size and impingement limit is required. An appreciation 
for the sensitivity, as well as the judgment required in the selection of 
size of cloud measuring elements and their application, can be obtained 
from the following example. The following table presents a spectrum of 
droplet sizes and the corresponding values of parameters and dimensions 
that apply: 


[Flight speed, 200 mph; altitude, 10,000 ft; 
airfoil chord length, 18 in. ] 


Droplet 

Inertia 

Free -stream 

Surface 

distance of impingement, s 

diameter, 

parameter, 

Reynolds 

Ratio to 

In. 

In inches for droplet 


K 

number. 

chord 


diameter difference 

microns 


Re 0 

length 


of 2 microns 

50 

1.64 

251 

0.230 

4.15 

) 

0.07 

48 

1.51 

240 

.227 

4.08 

/ 

40 

1.05 

201 

.206 

3.71 


.09 

38 

.95 

190 

.201 

3.62 

j 

30 

.59 

150 

.175 

3.15 


.14 

28 

.52 

140 

.167 

3.01 

/ 

20 

.26 

100 

.128 

2.30 

) 

.23 

18 

.21 

90 

.115 

2.07 

) 

10 

.066 

50 

.050 

.90 


.32 

8 

.042 

40 

.032 

.58 

I 

.17 

6 

.024 

30 

.023 

.41 

I 

j 

a .05 

4 

.010 

20 

.020 

.36 

1 

2 

.003 

10 

.018 

.32 

J 

a .04 


Obtained from extrapolated curves in fig. 8. 


The extent of impingement is given in the fourth column as a ratio 
of the actual distance to the chord length (obtained from fig. 8) and in 
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the fifth column in inches for an airfoil with an 18-inch chord. The 
last column gives the difference in limit in inches for a 2 -micron differ- 
ence in diameter of the droplets. In the table the largest difference in 
the limits and, therefore, the largest sensitivity, occurs between 10- 
and 8 -micron droplets although good sensitivity exists from 8- to 25- 
micron diameters. 

In the discussion of the dye-tracer technique presented in reference 
2, l/8-inch segments of blotter are used in the colorimetric analysis. 

Thus, the sensitivity is limited to variations that take place on incre- 
ments larger than l/8 inch. In the example cited in the preceding table, 
6-micron droplets cannot be differentiated from 4-micron droplets, and 
38-micron droplets cannot be differentiated from 40 -micron droplets when 
l/8-inch segments are used. A size spectrum would be most accurate between 
38 and 6 microns. If l/l6-inch segments are used (see ref. 10) , some im- 
provement is realized in the lower size range of the spectrum, but the 
lower limit is not extended appreciably , As can be seen from the preced- 
ing table, if the surface graduations were l/l6 inch apart, the lower 
range sensitivity would be extended to include the 5 -micron droplets. 

This extent ion would be an improvement of 1 micron in diameter. However, 
the use of l/l6-inch segments, especially near the stagnation line, im- 
proves the accuracy in determining the local rate of water distribution 
(Wp curve described in appendix B) . 

The limit of accuracy and sensitivity in the size spectrum (e.g. , 

38- to 6 -micron range in preceding example) will change when operating 
and meteorological conditions change. As discussed in the preceding sec- 
tion, a change in airfoil chord length has a large effect on the limits 
of sensitivity. The small airfoils 'have a desired advantage in measuring 
small droplets; however, decreasing the airfoil size introduces other 
difficulties. Decreasing the chord length from 18 to 6 inches decreases 
the surface distance (scale) on which the measurements are made, with a 
consequent requirement of greater care and more sensitivity. 

The preceding principles, actions, and conclusions apply equally well 
to cylinders and other geometric shapes. It should be noted that a cylin- 
der has different limits of sensitivity than an airfoil. Also, there are 
differences in sensitivity among airfoils of varying thickness ratio and 
shape as well as chord length mentioned in the preceding paragraph. 

Masking effect of large droplets . - In the determination of droplet - 
size distribution and liquid-water content, the masking effect of large 
droplets is a source of error inherent with most instruments that depend 
on impingement and retention of the cloud substance such as dye or ice 
accretion. With these instruments the final measurement, for example, 
the local rate of impingement on the surface of the element or total 
collection efficiency, is a composite of contributions by droplets of 
different sizes. The contribution from large drops may overwhelm the 
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contribution from the small droplets even though a large percentage of 
the total water in the cloud may he carried by the small droplets , because 
of the low collection efficiency of small droplets. This masking effect 
is illustrated with examples in appendix C obtained from table II. 

In appendix C two widely differing droplet-size distributions 
(flights 1 and 3, table II) measured with an aeroscope (ref. 3) are taken 
as the standards with which comparisons are made. These distributions 
are applied to a 36 .5-percent-thick Joukowski airfoil section used in 
such a manner that the local collection rate on chordwise sections of the 
airfoil surface is the principal measurable quantity. This method is used 
in the dye -tracer technique described in reference 2 and reviewed in ap- 
pendix B herein. The droplet-size distributions are divided into size 
groups. The water contained in each size group (fractional part of total, 
as measured with aeroscope) is given in rows D and K for flights 1 and 3, 
respectively. 

The. example evaluations for flight 1 are made at five chordwise 
positions on the airfoil surface (approximately 1 in. apart (rows E to I)) 
in order to present a condensed table of calculations . (The curve for the 
local rate of water impingement as a function of surface distance (see 
appendix B) is established by a large number of analyses of punched seg- 
ments (refs . 2 and 10) . For flight 3 the evaluations are made at the 
stagnation line only. 

Row E is an analysis at the airfoil stagnation line (s = 0) for 
flight 1. The rate of local droplet impingement p obtained from fig- 
ures 12 and 13 is given along with the rate of water impingement on the 
airfoil contributed by each size group (product of p • and water contained 
in each size group (row D) ) . The sum of the rates of water impingement 
contributed by each size group is the weighted local impingement effi- 
ciency. This value of weighted local impingement efficiency is 0.5482 at 
the stagnation line for flight 1. When the weighted value of p is used 
in equation (ll), the local rate of water impingement Wp is obtained. 

This value is used for the Wp curve in the dye -tracer technique for 
determining the droplet-size distribution in the cloud (see appendix B) . 

Errors in the measurement of dye concentration in each punched seg- 
ment will change the Wp curve used in determining the droplet-size dis- 
tribution in a proportional manner. A change in the Wp curve can affect 
the reported droplet-size distribution appreciably among the small drop- 
lets. Because the error in measuring dye concentration shrinks the Wp 
ordinate and not the surface extent of impingement and because the droplet- 
size distribution obtained by the dye-tracer technique is a cumulative 
composite of contributions by different sizes starting with the largest 
size, the small droplets contributing a small fraction of the total dye 
can be indiscernible in the analysis. For example, an error of -5 percent 
in measuring the dye concentration lowers the Wp by 5 percent. Through 
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equation (ll) the lowering of Wp can be considered in terms of a 
weighted local impingement efficiency p lowered by 5 percent. At the 
stagnation line in the example given in appendix C for flight 1 (row E) , 

5 percent lowers the value of weighted p by 0.0274. All the water 
impinging on the blotter at the stagnation line can be accounted for in 
the droplets between diameters of 57.5 and 12.6 microns. The amount of 
dye lost by the error is the same as that which is contributed by the 0- 
to 12. 6 -micron-diameter droplets. As can be seen by the analyses made in 
rows F to I, other positions on the airfoil surface do not aid in discern- 
ing these small droplets. When the droplets below 12.5 microns are in- 
discernible, the liquid-water content is only 0.89 as large as that 
measured by the aeroscope. 

The analysis along the stagnation line for flight 3 (row L) shows 
that a 5 -percent error obliterates the droplets between 0 and 7.5 microns. 
Elimination of these droplets makes the liquid-water content 0.95 as large 
as that measured with the aeroscope. 

A safeguard against losing too many of the small droplet groups is 
provided in reference 10 by using an aspirator to determine liquid-water 
content and, also, to establish the local rate of water impingement curve 
(appendix B) at the airfoil stagnation line. ' Furthermore, the use of 
several airfoils of different sizes (ref. 10) improves the inherent sen- 
sitivity on local rates of impingement as was shown for limits of 
impingement . 


Cylinder Surface Velocities 

The problem of discrepancies in cloud measurements when the size of 
the cylinder was varied is mentioned in the IMTR0DUCTI0W . According to 
the trajectory calculations in reference 1 based on the classical theory 
of airflow around a right circular cylinder, final cloud measurements 
(droplet-size distribution and liquid-water content) should not be depend- 
ent on the size of the cylinder used. The cylinder-size trend experienced 
in the work reported in reference 2, in which the final results were 
dependent on the particular cylinder diameter used, is caused by factors 
other than the sensitivity discussed in the preceding sections. Further- 
more, the trend is not part of the concept presented in reference 1 in 
which a group of cylinders of different sizes can be used for making 
cloud measurements . 

Because of the different surface to free-stream velocity ratio meas- 
ured on different-sized cylinders in the icing wind tunnel (fig. l), it 
was suggested in reference 2 that the differences in cloud measurements 
obtained by different sizes of cylinders were due to differences in the 
flow fields surrounding the cylinders . If the point of view is taken 
that the difference in surface velocities between the 6- and 2-inch 
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cylinders (fig, l) is caused "by the degree of separation of flow from the 
cylinder surface (fig. 16) and consequently by the location on the sur- 
face of the two cylinders where separation begins, the difference in im- 
pingement among the cylinders can be explained in part on the basis of 
differences in the flow field. When flow separation occurs, the airflow 
ahead of the cylinder is not that attendant to a cylinder but rather to 
some virtual shape such as shown in figure 16 by the forward part of the 
cylinder and the dashed lines . Since the trajectories of droplets are 
governed by the environmental airflow pattern, the impingement is related 
to the virtual shape caused by the flow separation. 

Icing tunnel experiments beyond those reported in reference 2 have 
shown that size trends were nearly eliminated when the cylinders were 
replaced by 36 .5 -percent -thick Joukowski airfoils of different sizes. 

The elimination of the size trend was due to the elimination of flow 
separation as can be surmised from the surface velocity distribution 
shown in figure 2. 

The differences in curves of surface velocity shown in figure 1 for 
different -sized cylinders may not apply under all operating conditions. 

A safe procedure in the use of cylinders is to survey the surface velocity 
on the cylinders with pressure taps . In order 'for a reasonable assurance 
that the theoretical trajectory results of reference 1 apply, the velocity 
survey should not reveal flow separation such as is revealed in figure 1. 

A pressure tap survey of surface velocities with the cylinders mounted on 
an airplane was made in flight in order to determine whether the separa- 
tion deduced from figure 1 was largely affected by tunnel conditions such 
as tunnel turbulence. In order to minimize end effects, both the 6- and 
2-inch flight cylinders were 15 inches long, and the holes for the pres- 
sure taps were located midway between the ends. 

The flight results are shown in figure 17 with the theoretical and 
6 -inch cylinder curves of figure 1 taken in the icing tunnel. The flight 
results are very comparable to the tunnel results in that the velocities 
of the 2-inch cylinder are lower than for the larger cylinder in both 
cases, and all are lower than the theoretical. The flight surface veloci- 
ties on the 4.5-inch cylinder of the multicylinder set detailed in figure 
16 of reference 1 are also shown in figure 17. Although the maximum ve- 
locity measured on the 4.5-inch cylinder is slightly higher than on either 
the 6- or 2-inch cylinder, the velocity everywhere is considerably lower 
than the theoretical and, thus, lower than the values necessary for safe 
application of the theoretical computations of reference 1. The somewhat 
higher maximum velocities obtained on the 4.5-inch cylinder mounted as 
part of the multicylinder set are not readily explainable. However, the 
fact that this cylinder was only about 3 inches long as compared to 15 
inches for the 6- and 2-inch-diameter cylinders and that the mounting was 
slightly different may have a bearing on an explanation. The 4.5-inch 
multicylinder was mounted as part of the multicylinder set. 
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Because the flight and tunnel velocity distributions on the surface 
of the cylinders are similar, the droplet impingement on the flight cylin- 
ders can be expected to behave in a manner similar to that on the tunnel 
cylinders . 


CONCLUDING REMARKS 

Viscosity was one of the effects mentioned in the INTRODUCTION that 
are not accounted for in theoretically derived flow fields. Although at 
present there is no direct experimental evidence evaluating the effect of 
viscosity on trajectories, there are reasons to believe that with most 
well designed aerodynamic bodies the effect of viscosity is small. An 
analytical evaluation of viscosity presented in reference 11 supports this 
point of view for the range of Reynolds numbers usually encountered in 
making natural cloud measurements from aircraft. Perhaps a simplified 
physical explanation of the reason the effect of viscosity on droplet 
trajectories around cylinders is negligible may be parallel to the expla- 
nation on why compressibility has small effect. The velocity field about 
the cylinder is influenced only near the cylinder surface by either com- 
pressibility or viscosity. Reference 1 shows for compressibility that 
when the flow field is affected only near the body, the effect is small on 
trajectories of droplets . 


The replacement of cylinders by a Joukowski airfoil for measuring 
instruments is required only when the airflow around the cylinder is not 
comparable to the theoretical flow used in the calculations of reference 1. 
The 36 .5-percent-thick Joukowski airfoil retains many of the desirable 
properties of the cylinder and has the necessary flow field that permits 
the application of calculated results in a range of operating conditions 
not possible with cylinders. The least desirable feature of the airfoil 
is that the shape changes as ice accumulates on the surface. This feature 
is of no consequence when the airfoil is used with the dye-tracer technique 
in tunnels . 

As mentioned in the INTRODUCTION, the impingement data presented 
herein are also useful in estimating impingement on thick streamlined 
struts. Thick-sectioned struts that are over 15 percent thick probably 
do not differ from Joukouski airfoil sections in any manner that affects 
the collection efficiency appreciably when compared to a Joukowski airfoil 
section of the same thickness ratio. The surface extent of impingement is 
more sensitive than collection efficiency to details of shape (ref. 10) . 

As demonstrated in reference 12, if a part of an airfoil section (such as 
the leading-edge region or forward 25 -percent region, etc.) has the same 
shape and pressure distribution as its counterpart on another airfoil, the 
impingement on the two parts will be the same even though the other por- 
tions of the airfoil sections differ in both shape and pressure 
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distribution. On thick-sectioned bodies the distribution is concentrated 
on the leading edge, which eases the problem of estimating the surface 
impingement . . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 9, 1957 
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APPENDIX A 
SYMBOLS 

a droplet radius, ft (3.048><10 5 microns) 

Cp coefficient of drag, dimensionless 

d droplet diameter, microns (3. 28xi0“® ft) 

E collection efficiency, dimensionless 

K inertia parameter, 1.704x10"-*-^ dimensionless (density of water, 

pL 

1.94 slugs/cu ft, included in constant) 

Kq modified inertia parameter, JL K, dimensionless 


L 

Re 


Re 0 


s 


t 

U 

u 

V 

v 

W 

W P 

w 

x,y 


airfoil chord length, ft 

local Reynolds number with respect to droplet, dimensionless 


free-stream Reynolds number with respect to droplet, 
dp U 

— dimensionless 
M- 


4. 813X10 “ 6 


distance on surface of airfoil measured from leading-edge chord 
point, ratio to chord length 

time, sec 

flight speed, mph 

local air velocity, ratio to free-stream velocity 
free-stream velocity, ft/sec 

local droplet velocity, ratio to free-stream velocity 

rate of water impingement per unit span of airfoil, lb/ (hr) (ft span) 

local rate of water impingement, lb/ (hr) (sq ft) 

liquid-water content in cloud, g/cu m 

rectangular coordinates, ratio to chord length 


QTrVfc 
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3 local impingement efficiency, dimensionless 

9 central angle of cylinder measured from leading edge stagnation 
point, deg 

X true range of droplet as a projectile injected into still air, ft 

X s range of droplet as a projectile following Stokes ' law, ft 

U viscosity of air, slugs/ (ft) (sec) 

p density, slugs/cu ft 

T dimensionless time functions, % = tV/L 

Reo 18p|LV 

<P = —rr- , dimensionless 

K UP W 

Subscripts : 
a air 

l lower airfoil surface 

max maximum or total rate 

s airfoil surface 

tan tangent trajectory 

u upper airfoil surface 

w water 

x,y component coordinate direction 

3 local rate 


0 


free stream 
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APPENDIX B 


REVIEW OF DYE-TRACER TECHNIQUE 

The method is reviewed herein to present the "background necessary 
for evaluating the use of the 36. 5 -percent -thick Joukowski airfoil as a 
sensing element in a cloud droplet instrument. The dye -tracer technique 
was developed in order to measure the droplet-size distribution and 
liquid-water content in a cloud. A "blotter-wrapped body , such as a cylin- 
der or the Joukowski airfoil studied herein, is exposed to an airstream 
containing a dyed-water spray cloud. In this technique, water treated 
with small quantities of water-soluble dye is injected in the form of 
droplets into the airstream ahead of the body by means of spray nozzles . 

At the point of impact and droplet absorption by the blotter, a permanent 
dye deposit is obtained. The amount of dye deposit on the blotter is 
directly proportional to the water impingement rate. The amount of dye 
trace obtained in a measured time interval can be determined by colorimet- 
ric analysis and converted into the quantity of water which produced it 
because the composition of the treated water is known. The local rate of 
water impingement Wg, total rate of water interception W max , and the 
surface extent of impingement s max are obtained from the colorimetric 
analysis. 


The principal information obtained from the colorimetric analysis is 
the local rate of water impingement W^ as a function of surface distance 
s. Typical values of W^ are shown in figure 18 as the solid-line curve. 
The area under the curve is a measure of the total water impingement on the 
blotter. The area under the curve is a composite of areas contributed 

by an unknown number of different droplet sizes in the cloud-size distribu- 
tion. The object is to determine the cloud-size distribution from the 
shape of the solid-line curve in figure 18. The determination is made with 
the use of the calculated 6 curves (local impingement efficiency) shown 
in figures 12 and 13. 

For simplicity in the presentation of the principles involved, the 
abscissa in figure 18 (airfoil surface distance of impingement) is arbi- 
trarily divided into six equal increments. As is described in reference 
2, the number and size of increments are usually governed by the size of 
the segments punched from the blotter for use in the colorimetric analysis 
(see fig. 12, ref, 2) . Hie limit of impingement s max is established by 
the largest droplets present in the droplet-size distribution. The actual 
diameter can be obtained from the values of s max and impingement limit 

data presented in figures 8 and 9, The contribution of these largest 
droplets to the total area under the solid line is obtained by superposing 
in figure 18 the proper (3 curve (local rate of droplet impingement) ob- 
tained from figures 12 and 13. This superposition properly weighted for 
a fractional part of total water is shown in figure 18 as the dashed line 


aTfYf? 
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labeled 1. If the solid-line curve were a result of impingement by a 
cloud composed of droplets uniform in size, curve 1 would coincide with 
the solid curve. However, figure 18 indicates that the largest droplets 
in the cloud contribute only the area under curve 1. 

The next step is to determine the contribution by a smaller size 
group. This group is assumed to have a limit of impingement at s 2 . 
Again, from the value of s 2 and the data of figures 8, 9, 12, and 13, 
the area contribution by group 2 is found to be the area between the 
dashed lines labeled 2 and 1. The procedure is repeated for S3, and so 
forth. Detailed procedure with working charts is given in reference 2. 
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APPENDIX G 


EVALUATION OP MASKING EFFECT 

An evaluation of the errors that actually are involved in the masking 
effect can he made with a realistic cloud distribution. The cloud data 
used in the following evaluation were obtained in flights through atmos- 
pheric clouds and are reported in reference 3. The measurements were 
made with an oil-stream photomicrographic aeroscope in which the cloud 
droplets were continuously captured in a stream of oil and photographed 
by a photomicrographic camera. 

The basic data used in table II (number of droplets in each size 
group) are obtained from the tables of reference 3. Two samples of size 
distributions taken on different days and in different cloud formations 
are used in table II and are identified with the flight, run, and picture 
numbers given in reference 3. The droplets are distributed by size 
(diameter) into groups, each group except the first covering a size range 
of 5 microns . The inertia parameter K (row A) and free-stream Reynolds 
number ReQ (row B) are calculated using the average diameter in the 
size group, except the first group for which -a 5 -micron diameter was 
used. 


The number of droplets in each size group is given in rows C and J 
for flights 1 and 3, respectively, and the amount of water in each size 
group, as a fractional part of the total water measured in the cloud, is 
given in rows D and K. For example, for flight 1 slightly more than one- 
tenth of the total water measured in the cloud is contained in the 269 
droplets in the 7 . 6 -to-12 . 5 -micron-diameter group. For flight 1 the rate 
of droplet collection is analyzed at five chordwise positions approximately 
1 inch apart on the surface of the airfoil (rows E to l) . For flight 3 
the analysis is made at the stagnation line only (s = 0, row L) . Each 
row gives the local collection efficiency p obtained from figures 12 
and 13 and the product of 3 and the water contained in the size group 
(rows D and K for flights 1 and 3, respectively), For example, at the 
stagnation line (s = 0) for flight 1 the local collection efficiency is 
0.19 for the 7.6- to 12. 5 -micron size group, and the product of p and 
the water contained in that size group is 0.19 x 0.107 = 0.0203. This value 
is the contribution by the 7.6- to 12. 5 -mi cron size group to a small seg- 
ment of blotter located at the stagnation line. The summation of all the 
contributions at the stagnation line by all the size groups present in 
the cloud gives a local collection efficiency which is weighted according 
to the droplet size versus the water content distribution in the cloud. 

The local rate of water impingement is obtained when the weighted value 
of p is used in equation (ll) . 
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In the dye-tracer technique the problem is worked, in the reverse 
order from that presented herein. The local rate of water impingement 
W p is the quantity measured, and from Wp the droplet-size distribution 
is obtained (see appendix B) . The examples in table II illustrate the 
sensitivity required in the dye-tracer technique in order to treat ade- 
quately size groups that make small contributions to the total effect. 

An error of -5 percent in p caused by errors in measurements of Wp 
is listed in each of the rows analyzed. The manner in which this error 
affects the droplet -size distribution and the liquid-water content is 
discussed in the text. 

For the purpose of the example illustrating the masking effect, it 
is assumed that the data on size distribution and liquid-water content 
measured by the aeroscope are as accurate as needed here, and the values 
measured by the aeroscope are arbitrarily taken as the standards for com- 
parison with an instrument which uses a Joukowski airfoil section (36.5 
percent thick) as the collecting element. 
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TABES I. - COORDINATE POINTS FOR SYMMETRICAL 
36. 5 -PERCENT JOUKOWSKI AIRFOIL 


X 

y 

s 

X 

y 

s 

0 

0 

0 

0.2803 

0.1818 

0.368 

.0024 

.0248 

.025 

.2884 

.1817 

.376 

,0095 

.0492 

.050 

.3713 

.1752 

.458 

.021 

.0726 

.078 

.4563 

.1594 

.548 

.0372 

.0945 

.104 

.5408 

.1369 

.633 

.0575 

.1146 

.132 

.6229 

.1106 

.723 

.0816 

.1325 

.161 

.7008 

.0834 

.800 

.1094 

.1479 

.194 

.7729 

.0580 

.879 

.1404 

.1606 

.227 

.8376 

.0364 

.948 

.1741 

.1704 

.262 

.8934 

.0199 

1.003 

.1884 

.1734 

.275 

.9388 

.0088 

1.049 

.2029 

.1760 

.291 

.9724 

.0027 

1.082 

.2104 

.1771 

.299 

.9931 

.0004 

1.103 

.2400 

.1810 

.328 

1.0000 

o 

1.110 




TABLE II. - EXAMPLE ILLUSTRATING MASKING EFFECT 
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Local surface to free- stream velocity ratio 







Velocity on airfoil surface, U B , ratio to free-stream velocity 
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Surface distance, s, ratio to chord length 


Figure 2, - Comparison of experimental (icing tunnel) with theoretical surface to free-stream 
velocity ratio for 36,5-percent-thick symmetrical Joukowski airfoil. Angle of attack, 0°; tun- 
nel airspeed, 175 miles per hour; airfoil chord length, 16,32 inches. 
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Figure 5* - Relation between distance along surface of 36.5-percent-thick Joukowskl airfoil 
and chordwise station. 



Inertia parameter 



(a) Free- stream Reynolds number , 16 . 

Figure 6. - Trajectory starting ordinates as function of point of impingement on surface of 36* 5-percent- thick Joukowski 
airfoil. Angle of attack, 0°. 



NACA TN 4035 


33 



'(b) Free- stream Reynolds number , 64. 

Figure 6. - Continued* Trajectory starting ordinates as function of point of impingement on surface of 36 . 5- percent -thick 
Joukovski airfoil. Angle of attack, 0 . 




Inertia parameter 



Figure 6. - Continued. Trajectory starting ordinates as function of point of impingement on surface of 36 . 5-percent- thick 
Joukovski airfoil. Angle of attack, 0 . 







Inertia parameter 



Figure 6. - Concluded. Trajectory starting ordinates as function of point of impingement on surface of 36. 5-percent- thick 
Joukowski airfoil. Angle of attack * 0 . 
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0 .004 .008 .012 .016 .020 .024 .028 

Distance on airfoil surface, s, ratio to chord length 


(c) Inertia parameter, 0.02008. 

Figure 7. - Trajectory starting ordinates as function of point of impingement on airfoil surface for 
small values of inertia parameter . Angle of attack, 0°} 36 .5 -percent -thick Joukovski airfoil. 






Limit of impingement on airfoil surface, a mA Y> ratio to chord length 
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Inertia parameter , K 
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Figure 11. - Collection efficiency of 36.5-percent-thick Joukovski airfoil as function of inertia 
parameter and (p. Angle of attack, 0°. 
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(a) Free-stream Reynolds number, 16. 

Figure 12. - Local impingement efficiency on surface of 36 .5-percenfc-thick Joukowski airfoil. Angle of 
attack, 0°. 
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(b) Free-stream Reynolds number, 64. 

Figure 12. - Continued. Local impingement efficiency on surface of 36 . 5-percent-thick Joukowskl air- 
foil. Angle of attack, 0°. 
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(e) Free -stream Reynolds number, 128. 


Figure 12. - Continued. Local Impingement efficiency on surface of 36.5-percent-thick Joukowski air- 
foil. Angle of attack, 0°. 
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0 .04 ,08 .12 .16 .20 .24 .28 .32 .36 .40 

Distance on airfoil surface, a, ratio to chord length 

(e) Free-stream Reynolds number, 1024. 


Figure 12. - Concluded. Local impingement efficiency on surface of 36.5-percent-thick Joukowski air- 
flow. Angle of attack, 0°. 
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(a) Inertia parameter, 0.08035. 
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Figure 14. - Comparison of collection efficiency for 15- and 36.5-percent-thick Joukovski airfoils with 
cylinder . 
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Figure 15 ^ - Comparison of surface extent of impingement for 15- and 
36.5-percent-thick: Joukowski airfoils with cylinder. 
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Airfoil surface distance, s 

Figure 18. - Illustration of total impingement rate on airfoil surface with con- 
tributions from several droplet size groups. 
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